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物 特异 性 ， 利 用 qRT-PCR 技术 检测 各 个 
GeNorm, NormFinder, BestKeeper 软件 和 RefFinder 在 
TFL1 的 表达 分 析 验 证 其 可 靠 性 。 结 果 表 
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研究 中 稳定 表达 的 内 参 基因 。 该 研究 结果 将 为 木兰 属 植 物 二 次 成 花 分 子 调控 机 制 研究 提供 
依据 。 
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Abstract: A number of Magnolia (Magnoliaceae) species bloom twice each year, instead of 
once in most other species in this family, which is a desirable ornamental trait. To investigate 
the molecular mechanism of the flower bud differentiation during the second bloom each year in 
these Magnolia species, quantitative real-time PCR (qRT-PCR) was frequently used as a 
sensitive gene expression technique that relies on the stability of reference genes for data 
normalization. Therefore, the identification of reference gene(s) suitable for molecular 
characterization during Magnolia flower bud differentiation in second bloom is of great interest. 
In this study, gene expression in the leaf and flower bud tissues of M. liliflora *Hongyuanbao' at 
different stages of the flower bud differentiation during second bloom was analyzed. Based on 
transcriptomic sequencing data, eight constitutively expressed genes, including UBC (ubiquitin- 
conjugating enzyme), ACT (actin), -TUB (tubulin beta), -TUBS (tubulin beta), a-TUB3 
(tubulin alpha), PEPC (phosphoenolpyruvate carboxylase), 4CP2 (acyl carrier protein 2), ACP3 
(acyl carrier protein 3), were selected as candidate reference genes for qRT-PCR. 
Comprehensively analysis was conducted using four softwares including geNorm, NormFinder, 
BestKeeper and RefFinder. Primer Premier 5 was used to design the primers. PCR products of all 
the eight candidate reference genes were analyzed by gel electrophoresis which showed sharp 
bands with the expected size. Each melting curve showed a single peak, which indicated the high 
specificity of PCR primers. The general assessment by the four different softwares ranked f-TUB, 
p-TUB5 and a-TUB3 as the most stable reference genes, whereas UBC and ACT were the lest 
stable. The reference genes were further evaluated by analyzing the relative expression of TFL/ 
gene with either single or in combination of £-TUB5, a-TUB3, B-TUB, which showed highly 
consistent results. In contrast, ACT and UBC did not effectively standardize the expression level 
of TFLI. In general, through this study, we have identified f-TUB5, a-TUB3 and f-TUB as the 
most suitable reference genes for qRT-PCR analysis of gene expression in flower bud 
differentiation during the annual second bloom of M. liliflora 'Hongyuanbao', which provides 


useful tools for investigating the regulatory mechanisms in Magnolia flowering. 
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木兰 属 CMagnolia Linn.) 植物 在 我 国 已 有 两 千 多 年 的 栽培 历史 ， 花 色 艳 丽 、 树 形 优 
美 ， 为 重要 的 春季 开花 绿化 植物 ( 张 缀 等 ，2009) 。 在 自然 条 件 下 ， 木 兰 属 植物 中 少数 
质 能 在 夏季 二 次 开花 ， 甚 至 可 持续 开放 到 秋季 (Xia et al., 2008) 。 然 而 ， 对 木兰 属 植物 
二 次 开花 或 多 次 开花 机 理 研 究 仍 处 于 初级 阶段 ， 目 前 仅 有 紫 玉 兰 “ 红 元 宝 ” (Magnolia 
liliflora *Hongyuanbao') (fé/b BÆ, 20200 和 常春 二 乔 玉兰 (Magnolia soulangeana 
*Changchun' ) (Sun et al., 20210 两 个 品种 不 同 花 期 的 转录 组 分 析 研 究 ， 从 基因 调控 层面 
揭示 该 属 植物 二 次 开花 机 制 的 内 容 鲜 有 报道 。 

实时 荧光 定量 PCR (quantitative Real-time polymerase chain reaction; qRT-PCR) 技术 
具有 高 精度 、 高 通 量 、 高 灵敏 度 等 优点 ， 被 广泛 用 于 基因 表达 分 析 的 研究 中 CHuggett et 
al., 2005) 。 但 qRT-PCR 结果 的 准确 性 容易 受到 样本 差异 性 、 组 织 特异 性 、 提 取 RNA Ji 
量 等 因素 的 影响 (Suzuki et aL, 2000) ， 需 要 引入 一 个 或 几 个 内 参 基因 校正 目的 基因 的 表 


Ek 


ik (Quackenbush, 2002; Xu et al., 2017) 。 实 


范围 内 表达 恒定 ， 不 能 确保 在 所 有 试验 条 件 
经 筛选 而 盲目 使 用 他 人 试验 中 的 内 参 基 因 ， 

al., 2006) 。 因 此 ， 在 不 
达 的 内 参 基 因 。 有 目前 
处 理 〈 王 倩 颖 等 ，2019) ZEB AE 


际 研究 中 所 选取 的 内 参 基因 通常 只 能 在 一 定 


样本 或 不 同 试验 条 


zi 


均 能 稳定 表达 CThellin et al., 1999) 。 若 不 


容易 得 到 不 准确 甚至 错误 的 结果 (Mukesh et 


件 下 进行 qRT-PCR 4 ITIN BU $e ze Ur E e XE 


已 有 研究 筛选 得 到 木兰 属 


开花 的 木兰 属 植物 进行 内 参 基因 的 筛选 。 

紫 玉 兰 “ 红 元 宝 ” 
(Magnolia liliflora) 的 变种 ， 具 有 一 年 二 次 开花 的 特性 〈 程 少 怠 等 ，2020) 。 为 筛选 适 
合 木 兰 属 植物 二 次 花芽 分 化 研究 的 内 参 基 因 ， 本 研究 以 “ 红 元 宝 ' 不 同 花 芽 分 化 时 期 的 花 


芽 、 叶 为 材料 ， 基 于 
a-TUB3, PEPC, ACP2, ACP3) ， 并 利 月 


C Magnolia liliflora ‘H 


植物 不 同 花 色 〈 王 宁 杭 等 ，2019) . Ra 


〈 常 鹏 杰 等 ，2018) 适用 的 内 参 基 因 ， 却 未 对 二 次 


[o 


ongyuanbao') (简称 “ 红 元 宝 ”) 是 紫 玉 兰 


转录 组 数据 选取 8 个 常用 内 参 基 因 (UBC、ACT、p-TUB、p-TUB5、 
H geNorm, NormFinder 和 Bestkeeper 等 统计 学 


软件 综合 评估 候选 内 参 基因 在 不 同 花 芽 分 化 时 期 的 表达 稳定 性 ， 以 期 为 木兰 属 植物 二 次 或 


多 次 花 开 放 的 分 子 机 到 


1 材料 与 方法 
1.1 材料 选取 


民 好 无 病虫害 的 “ 红 


“ 红 元 宝 ” 《树龄 


ER 


— ct 


-LJUE 


完 提 供 依据 。 


& 10 年 ) 种 植 于 浙江 农 


林 大 学 平山 基地 。 选 取 三 棵 长 势 一 致 ， 状 况 


为 试验 材料 ， 定 期 采集 树冠 中 上 部 分 的 花芽 及 其 周围 的 嫩 叶 。 


由 于 每 棵 树 上 花芽 分 化 的 进程 并 不 同步 ， 因 此 不 同 分 化 时 期 的 花芽 按照 对 应 的 尺寸 规格 进 
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j 12 日 进行 第 一 次 采样 ， 


每 隔 15 天 采样 一 次 ， 采 集 三 次 作为 第 一 次 


花芽 分 化 前 -后 三 个 时 期 的 样品 2019 年 6 月 25 日 开始 采集 第 二 次 花芽 分 化 样品 ， 每 隔 


15 天 采集 
MES 3^8 


1.2 5| 


根据 “ 红 元 宝 ” 转 录 组 
(ACT) 、 
求 ， 运 月 


物 设计 


BEL 2. 


13 实 


验方 法 


ER, H36 个 样品 。 
中 超低温 保存 ， 最 终 转 到 -80 'C 冰 箱 储 存 备 用 。 


的 测序 结果 ， 从 中 选取 泛 素 酶 基因 (UBC) . Wiz) A 
B 微 管 蛋 白 〈p-7UB ) 等 8 个 常见 的 内 参 基 因 。 根 据 菊 光 定 量 引 物 设计 的 要 
H Primer Premier 5 软件 设计 引物 。 引 


次 ， 采 集 三 次 作为 第 二 次 花芽 分 化 前 中 后 三 个 时 期 的 样品 。 每 个 时 期 的 样品 
每 个 样品 均 置 于 5 ml 的 冻 存 管 中 ， 采 集 后 迅速 放 入 液 氮 


物 交 于 浙江 有 康生 物 科技 有 限 公 司 ， 序 列 信 


按 UltraClean Polysaccharide and Phenol Plant RNA Purification Kit (天津 ，Novogene) 


试剂 盒 
整 性 、 


说 明 书 提取 各 检 


浓度 和 纯度 。 有 月 


进行 10 uL 体系 反 转 录 ， 


fF 品 的 总 RNA。 通 过 1% 琼 脂 糖 凝 胶 电 泳 和 紫外 分 光 光 度 计 检测 其 完 
肯 反 转录 试剂 盒 PrimeScriptTM RT Mix (TaKaRa) 进行 cDNA 合成 ， 
取 2 uL 5xPrimeScript 缓冲 液 ，500 ng 总 RNA, RNase Free dH2O 


补足 至 10 aL, T 37 CC 反应 15 min Ji, 85 CE 反应 5Ss， 最 后 在 4 保存 。 用 LightCycler 


480 II (Roche) X J6 7 i 
uL， 上 下 游 引物 各 0.8 uL, 25 ng cDNA 模板 ， 


次 技术 重复 。qRTPCR 扩 增 反应 程序 : 95'C 


EAX BEIT qRT-PCR 反应 ， 其 中 BCG qPCR Master Mix (2x) 10 


环 ， 后 运行 95 C 5s, 60 C 60s, 95 C 15s. 


1.4 数 


本 研究 基于 “ 红 元 宝 
(UBC) 、 肌 动 蛋白 CACT) 、 微 管 蛋 白 B 链 
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据 分 析 


”转录 组 数据 ， 选 ] 


首 用 双 蒸 水 补足 至 20 此 ， 每 个 样品 进行 3 
预 变 性 30 s，95 'C 5s, 60 C 30s，40 个 循 


取 了 8 个 常见 的 候选 内 参 基 因 : 泛 素 酶 基因 
(B-TUB) ~ WEE AB-5 SE (B-TUB5) 、 微 


ÁH a-3 链 (a-TUB3) 、 磷 酸 烯 醇 丙 本 


HR RMI (PEPC) 、 栈 基 载 体 蛋白 2 


CACP2) 、 栈 基 载 体 蛋 白 3 (4CP3) 。 通 过 ACt，geNorm，NormFinder，Bestkeeper， 
RefFinder 分 别 计算 qRT-PCR 中 所 得 到 的 数据 ， 对 其 稳定 性 进行 稳定 性 分 析 ， 最 终 确 定 最 
佳 的 单个 内 参 基因 或 内 参 基因 组 合 。 候 选 内 参 基因 的 表达 水 平 通常 是 由 qRT-PCR I0] B 4 
周期 Ct 值 决定 ，Ct 值 与 基因 表达 水 平成 反比 ，Ct 值 越 大 ， 基 因 表 达 水 平 越 低 ， 反 之 亦 
然 。geNorm 是 根据 基因 表达 稳定 值 CM 值 ) 的 大 小 来 判断 候选 内 参 基 因 的 稳定 性 
(M«1.5) , M 值 越 小 意味 着 基因 表达 越 稳 定 。 为 了 确定 能 在 一 定 范 围 内 恒定 表达 的 最 优 
内 参 基 因数 量 ，geNorm 则 能 通过 配对 差异 值 Casa. 来 判断 最 佳 内 参数 量 从 而 进行 归 
化 (Vandesompele et aL, 2002) 。 若 基因 间 成 对 差异 值 大 于 或 等 于 0.15， 则 应 多 增 一 个 基 
因 对 其 进行 可 靠 的 归 一 化 。 帮 成 对 差异 值 小 于 0.15， 则 不 需要 增加 基因 进行 正常 化 
(Liang et al., 2020) 。NormFinder 将 原始 Ct 值 进行 了 对 数 转换 ， 然 后 进行 方差 的 分 析 ， 
计算 得 出 一 个 基因 的 稳定 表达 值 ， 稳 定 值 最 小 的 候选 内 参 基 因 则 为 最 稳定 的 基因 
(Lindbjerg et al., 2004) > BestKeeper 通过 计算 Ct 值 的 变异 率 比 较 不 同 候选 内 参 基 因 的 表 
达 水 平 差异 ， 相 关系 数 也 越 大 ， 标 准 偏差 SD 和 变异 系数 CV 越 小 则 越 稳定 (Pfaffl et al., 
2004) 。RefFinder 是 一 个 用 于 评估 和 筛选 参考 基因 的 在 线 网 站 ， 它 集成 了 geNorm， 
NormFinder 和 BestKeeper 的 计算 程序 ， 按 照 每 个 程序 对 候选 内 参 基因 的 排名 ， 为 每 个 基 
因 分 配 适当 的 权重 ， 并 计算 出 它们 的 权重 几何 平均 值 ， 从 而 获得 最 终 的 总 体 排名 〈 黄 真 池 
等 ，2013; Xie etal., 2012) 。 
2 结果 与 分 析 
2.1 总 RNA 质量 检测 

总 RNA 的 琼脂 糖 凝 股 电泳 检测 所 得 条 带 清晰 ， 并 且 28S rRNA 条 带 亮度 约 18S rRNA 
的 2 倍 ， 无 严重 的 RNA 降解 和 DNA 污染 情况 (图 1) 。 紫 外 分 光 光 度 计 检 测 结果 显示 ， 
所 有 RNA 样品 中 OD 260/OD zg0 值 均 在 1.8 ~ 2.0 £A, WH RNA 纯度 较 高 ， 可 用 于 进 一 
步 研 究 。 
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M. Marker; 28S. 28S rRNA; 18S. 18S rRNA; 5S. 5S rRNA; F1-F6. 两 次 花芽 分 化 6 个 时 期 的 花 
Zi L1~L6. 两 次 花芽 分 化 6 个 时 期 的 叶 。 
M. Marker; 28S. 28S rRNA; 18S. 18S rRNA; 5S. 5S rRNA; F1 ~ F6. Flower buds at six stages of two flower 


bud differentiation; L1 ~ L6. Leaves at six stages of two flower bud differentiation. 
图 1 ATE’ 总 RNA 琼脂 糖 凝 胶 电 泳 检 测 

Fig.l Agarose gel analysis of total RNA extracted from Magnolia liliflora *Hongyuanbao' 
2.2 候选 内 参 基因 引物 评价 

将 “ 红 元 宝 ” 不 同 花 芽 分 化 时 期 样品 的 cDNA 等 量 混 合作 为 模板 ，8 个 候选 内 参 基因 
通过 简单 克隆 扩 增 之 后 ， 各 候选 内 参 基因 结果 均 与 预期 一 致 ， 条 带 均 显示 在 80 ~ 150 bp 
之 间 且 单一 清晰 (图 2) 。 随 之 将 cDNA 稀释 成 6 个 浓度 梯度 (55, 57, 52, 53, 54, 5- 
5) ， 以 不 同 浓度 的 cDNA 作为 qRT-PCR 扩 增 的 模板 。 通 过 绘制 标准 曲线 法 ， 以 不 同 浓度 
的 对 数值 为 横 坐 标 ，Ct 值 为 纵 坐 标 即 可 获得 斜率 (Kk) 和 相关 系数 (xr?) 。 运 用 扩 增 效率 


公式 : E= C101) x10095, n[AU EEIE 90.0% ~ 110.0% 之 间 ， 相 关系 数 r?20.98 (K 
T 


1) 。 各 候选 内 参 基因 的 单一 炊 解 峰 也 表明 引物 特异 性 良好 《图 3) ， 说 明 实验 结果 准 


可 靠 。 
2 000 bp 
1 000 bp 
750 bp 
500 bp 
250 bp 
100 bp 
M. Marker; 1. UBC; 2. ACT; 3. B- TUB; 4. B-TUB5; 5. a-TUB3; 6. PEPC; 1. ACP2; 8. ACP3. 
2 8 个 候选 内 参 基因 PCR 扩 增 产物 
Fig.2 PCR amplification products of 8 candidate reference genes 
表 1 8 个 候选 内 参 基因 的 引物 序列 及 扩 增 效率 
Table 1 Primers and amplification efficiency of 8 candidate reference genes 
基因 名 称 引物 序列 (5'—>3') 产物 长 度 bp ” 扩 增 效率 (%) 相关 系数 r 
Gene name Primer sequence Product length Efficiency Correlation coefficient 
UBC F_TAACAACGGCAAACCCAG 109 94.30 0.983 
R_GCAGTCGAATGCCTAATC 
ACT F_AGTCCAGCACAATACCAG 128 108.85 0.988 
R_CAAGGCTAATCGTGAGAA 
PB-TUB F_CCCGTCTCCACTTCTTCA 121 91.99 0.993 
R_CACATCATATTCTTAGCATCCC 
p-TUB5 F ATCAGTTCCGCTCCTTCA 81 103.80 0.995 
R. GCCCGATAACTTCGTCTT 
a-TUB3 F CTGGCTTGGGTTCTTTAC 100 97.80 0.993 
R. GCCGTTGATACCTGTGGG 
PEPC F GGGCATTGTAAGGGAGGC 121 101.09 0.996 
R. ATGCGGGCAGGAATGAGT 
ACP2 F GCAGCAGTGGGTAAGCAA 95 92.34 0.997 
R. GAGCGTGAACCAGAATAT 
ACP3 F CATAACAAGCTCCACCCT 141 99.08 0.995 


F ACCAGCTCCGATCATCTA 
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-dF/dT. Negative derivative ratio of fluorescent signal to temperature. 
图 3 8p VI RE DATA AE V 2 


Fig.3 Melting curve of 8 candidate reference genes 


2.3 内 参 表 达 稳 定性 分 析 
通过 比较 8 个 候选 内 参 基 因 Ct 值 来 评估 其 在 “ 红 元 宝 ”不 同 花 芽 分 化 阶段 的 表达 丰 
度 ， 根 据 Ct 值 的 分 布 (图 4) ，8 个 候选 内 参 基 因 的 Ce 值 均 集 中 在 17 ~ 30 之 间 ， 表 达 丰 


度 适 中 。 
24 
22 
20 


UBC ACT p-TUB f-TUB5 a-TUB3 PEPC ACP2 ACP3 
候选 内 参 基因 Candidate reference gene 
箱 体 表示 Ct 值 的 集中 范围 ， 箱 体 中 的 横 线 表示 中 位 数 ， 箱 体 上 下 边 分 别 表示 上 /下 四 分 位 数 ， 箱 体 上 下 
两 端 分 别 表 示 最 大 /小 值 ， 箱 体内 方 框 表示 均值 ， 箱 体外 黑 点 表示 变异 值 。 


The box represents the concentrated range of Ct values, the horizontal line in the box represents the median, the 


a 


Ct 值 Ct value 


X 


upper and lower sides of the box respectively represent the upper/lower quartiles, the upper and lower ends of the 
box respectively indicate the maximum/minimum value, the square box inside the box body represents the mean 
value, and the black dot outside the box body represents the variation value. 

图 4 8 个 候选 内 参 基因 的 Ct 值 


Fig.4 Ct value of 8 candidate reference genes 


geNorm 对 各 个 候选 内 参 基 因 进 行 稳定 性 排序 ， 结 果 显 示 (图 5) ，8 个 M 值 均 小 于 
1.5， 最 稳定 的 内 参 基 因为 B-TUB 和 p-7TUB5， 表 达 稳 定 值 最 小 (M=0.598 ) 。 UBC 


(M-1.137) 和 ACT (M=1.021) 在 该 软件 的 分 析 中 排序 靠 后 ， 是 最 不 稳定 的 内 参 基 因 。 
的 分 析 显 示 〈 图 6) , Vasa 的 值 均 大 于 0.15， 其 原因 可 能 是 “ 红 元 宝 ” 二 次 


配对 差异 值 
成 花 时 间 跨 度 大 ， 
个 默认 值 ， 并 非 总 


( 任 锐 等 


不 同 组 织 中 基因 变异 大 。 在 geNorm 手册 中 也 提 及 Van W 0.15 只 是 一 
是 必要 ， 可 根据 研究 目的 或 Vun 变化 趋势 选择 2 ~3 个 稳定 的 内 参 基 


，2016) . NormFinder 软件 分 析 基 因 稳 定性 同样 与 数值 六 小 呈 负 相关 ， a- 


TUB3 的 MM fi zk/ (C M=0.214) ， 因 此 稳定 性 最 高 ， 稳 定性 次 之 的 是 Bp-TUB5 
(M=0.308) ， 稳 定性 最 差 的 是 UBC (M=0.910) 和 4CT (M=0.847) ， 这 一 结果 与 
geNorm 大 致 相同 。BestKeeper 分 析 结 果 显 示 ACT 的 SD 值 最 大 (SD=2.15) ， 是 最 不 稳定 


的 内 参 基因 ， 


但 4CP2 的 SD 值 最 小 (SD=1.04) 、UBC 次 之 (SD=1.08) ,表明 ACP2 和 


UBC 的 稳定 


性 较 好 ， 这 与 geNorm、NormFinder 的 分 析 结 果 部 分 不 一 致 。RefFinder 集合 


了 上 述 三 种 软件 的 算法 和 程序 ， 对 每 个 备 选 基因 计算 稳定 值 重 的 几何 平均 值 进行 综合 排 
名 ， 数 值 越 小 ， 则 基因 越 稳定 〈 陈 莹 等 ，2019) 。 通 过 分 析 结 果 得 出 〈 表 2) , a-TUB3 
表达 最 为 稳定 ，p-7UB5 和 8-TUB 稳定 性 次 之 。 
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基因 表达 稳定 值 C(M 值 ) 
Stable value of gene expression (M value) 
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图 5 geNorm 软件 分 析 候 选 内 参 基 因 表 达 稳 定 值 


Fig.5 Stable value of gene expression of candidate reference gene analyzed by geNorm 
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图 6 geNorm 软件 分 析 候 选 内 参 基 因 的 Vasa 


Fig.6 Pairwise variation value (Vwn+1) of candidate reference gene analyzed by geNorm 


基于 上 述 3 个 软件 的 分 析 结 果 CgeNorm, NormFinder. Doyleeper? ， 通 过 RefFinder 
综合 确定 内 参 基因 的 稳定 性 排名 。 结 果 显 示 ( 表 2) ， 在 “ 红 元 宝 ” 不 同 花 芽 分 化 阶段 ， 
a-TUB3. B-TUB5. B-TUB 的 表达 稳定 性 最 高 ， 而 UBC 和 4CT 表 达 稳 定性 最 低 。 
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X2 候选 内 参 基 因 稳 定性 排序 


Table 2 Stability sequencing of candidate reference gene 


候选 内 参 基因 geNorm NormFinder BestKeeper RefFinder 综合 排序 
Candidate reference Comprehensive ranking 

UBC 1.137 0.910 1.08 5.66 V 

ACT 1.021 0.847 2.15 7.24 8 

p-TUB 0.598 0.496 1.55 3:13 3 
p-TUB5 0.598 0.308 1.61 2.30 2 
a-TUB3 0.718 0.214 1.46 1.97 1 
PEPC 0.910 0.509 1.36 4.95 6 

ACP2 0.843 0.562 1.04 3.66 5 

ACP3 0.781 0.451 1.15 3.22 4 


2.4 TFLI 基因 验证 最 佳 内 参 基 因 的 稳定 性 

TERMINAL FLOWER 1 (TFL1) 是 重要 的 成 花 抑 制 基因 ， 因 此 ， 选 用 TFL 基因 对 候 
选 内 参 基因 进行 验证 。 结 果 表 明 (图 7) , TELI 在 叶 中 几乎 没有 表达 ， 而 在 花芽 中 ， 第 
一 次 花芽 分 化 的 表达 量 明显 低 于 第 二 次 花芽 分 化 ， 推 测 低 水 平 表达 的 TELI 可 能 促进 了 
“ 红 元 宝 ” 二 次 开花 。 而 以 稳定 性 较 差 的 内 参 基因 分 析 TFLI 的 相对 表达 量 时 ，7ZZ7 基 
因 表达 模式 发 生 了 变化 。 以 单个 基因 (6-TUB 、p-TUB5、a-TUB3) 及 其 组 合 为 内 参 基因 
校准 TFLI 基因 的 表达 水 平时 ， 均 显示 出 了 一 致 的 趋势 ， 进 一 步 验证 了 4 个 软件 对 候选 内 
参 基 因 筛 选 结果 的 可 靠 性 。 


ni 
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第 一 次 花芽 分 化 前 期 Early stage of the first flower bud differentiation HESS 第 二 次 花 丁 分 化 前 期 Early stage of the second flower bud differentiation 
UA 第 一 次 花芽 分 化 中 期 Middle stage of the first flower bud differentiation E— 第 三 次 花芽 分 化 中 期 Middle stage of the second flower bud differentiation 
SSS 第 一 次 花芽 分 化 后 期 Late stage of the first flower bud differentiation HE 第 二 次 花芽 分 化 后 期 Late stage of the second flower bud differentiation 


A. 7TF7Z7 在 花芽 中 的 相对 表达 量 ，B. TFLI 在 叶 中 的 相对 表达 量 。 
A. Relative expression of TFL1 in flower buds; B. Relative expression of TFL/ in leaves; 1. UBC; 2. ACT; 
3. B-TUB; 4. B-TUB5; 5. a-TUB3; 6. PEPC; 7. ACP2; 8. ACP3; 9. B-TUB. p-TUB5. a-TUB3. 

7 不 同 内 参 基 因 或 组 合 分 析 TFL] 在 花芽 分 化 阶段 的 相对 表达 量 


Fig.7 Relative expression of 7FZ7in flower bud differentiation using different 
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reference genes or combinations 
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3 讨论 与 结论 


理想 的 内 参 基因 在 任何 试验 条 件 下 都 能 稳定 表达 ， 但 大 量 研 究 表明 ， 在 不 同 物 种 或 不 
同 条 件 中 稳定 表达 的 基因 会 发 生变 化 CArtico et al., 2010; Mafra et al., 2012) 。 因 此 ， 在 
特定 条 件 下 第 选 合适 的 内 参 基因 对 于 后 续 基 因 表达 模式 分 析 具 有 重要 意义 。 本 研究 通过 转 
录 组 测序 数据 ， 和 筛选 出 “ 红 元 宝 ” 花 芽 分 化 进程 8 个 候选 内 参 基因 ， 基 于 GeNorm、 
NormFinder、BestKeeper 和 RefFinder 软件 综合 分 析 结 果 ， 进 一 步 验 证 了 p-TUB、p-TUB5 
和 a-7TUB3 Æ ANR? 花芽 分 化 进程 最 佳 内 参 基 因 。 
E ALR? 不 同 花 芽 分 化 阶段 ，GeNorm、NormFinder 和 BestKeeper 软件 评估 的 8 
个 内 参 基 因 排 序 存 在 差异 ( 表 2) 。 有 研究 显示 ， 利 用 GeNorm. NormFinder 和 
BestKeeper 软件 分 析 莱 莉 花 《Jasminum sambac) (FEES, 20200 WIES (Amana 
edulis) (RIRE, 2021) 不 同 组 织 的 内 参 基因 表达 量 时 ，3 个 软件 得 到 的 稳定 性 排名 
也 不 同 。 由 此 ，GeNorm、NormFinder 和 BestKeeper 软件 分 析 结 果 的 差异 在 物种 中 普遍 存 
在 〈 付 建新 等 ，2016) 。 在 本 研究 中 ，GeNorm 和 NormFinder 分 析 结 果 表 明 a-TUB3、p- 
7TUB5 和 Bp-7TUB 稳定 性 较 好 ，UBC 和 4CT 表达 稳定 较 差 ， 而 BestKeeper 与 前 二 者 结果 不 
一 致 ， 认 为 ACP2 和 UBC 稳定 性 最 佳 ，4CT 稳 定性 最 差 ( 表 2) 。 造 成 分 析 结 果 有 差异 
的 原因 是 3 个 软件 采用 的 统计 学 算法 不 同 ， 使 同一 组 试验 数据 产生 了 不 同 的 排名 
(Spiegelaere et al., 2017) o 
为 了 减少 软件 算法 差异 造成 的 偏差 ， 通 常 结合 几何 平均 值 法 整合 不 同 软件 的 结果 得 出 
综合 排名 〈 任 锐 等 ，2016) 。 本 研究 利用 RefFinder 综合 评估 上 述 3 个 软件 的 排序 结果 ， 
分 配 权重 进行 几何 平均 值 的 计算 ， 由 表 2 可知， 在 所 有 样品 中 稳定 性 最 佳 的 是 6-TUB、p- 
TUBS5 和 a-TUB3， 稳 定性 最 差 的 是 UBC 和 A4CT。 为 进一步 验证 筛选 结果 的 可 信和 度 ， 本 而 
究 使 用 所 有 单个 内 参 及 最 佳 内 参 基因 组 合 (6-TUB、p-TUB5 flla-TUB3) 对 目的 基因 TFLI 
进行 相对 表达 量 分 析 CHE TO ， 结 果 表 明 3 个 内 参 基 因 组 合 校准 的 基因 表达 模式 更 为 准 
确 。 在 老 鸦 汶 (Amana edulis) GRBE, 20210 . KRIE (Jasminum sambac) (XF 
玉 等 ，2020) . © (Miscanthus sinensis) 〈 黄 逸 之 等 ，2020) 等 物种 的 内 参 研 究 中 也 证 
实 ， 综 合 不 同 软件 的 分 析 结 果 更 能 反映 给 定 试验 条 件 下 内 参 基因 的 表达 稳定 性 。 由 此 表 
= 明 ， 本 研究 利用 多 个 不 同 的 算法 对 基因 稳定 性 进行 综合 评估 从 而 对 候选 内 参 基 因 进 行 得 
- 选 ， 研 究 结果 准确 可 靠 。 
€) TUB 基因 是 生物 体 维持 生命 活动 所 必须 的 细胞 器 骨架 的 基本 组 分 ， 由 co-7UB 和 pf-TUB 
组 成 ， 具 有 和 良好 的 稳定 性 〈 阙 友 雄 等 ，2009) 。 有 研究 结果 表明 ， 马 继 杜 鹃 
(Rhododendron delavayi) 花芽 的 不 同 发 育 时 期 的 内 参 研究 中 ，7TUB 是 最 佳 的 内 参 基因 之 
一 (张明超 ，2019〉; WAE (Vitis davidii) CRAE, 2019) Mazra CLycoris) CR 
婷 婷 等 ，2015) 在 不 同 组 织 器 官 和 非 生物 胁迫 条 件 下 筛选 内 参 基因 时 ， 发 现 TUB 的 稳定 
性 最 好 。 与 前 人 研究 结果 一 致 ， 本 研究 综合 4 个 软件 分 析 结 果 确 定 TUB 同 源 基 因 p- 
7UB、p-7UB5 flla-TUB3 为 “ 红 元 宝 ” 二 次 花芽 分 化 研究 的 最 佳 内 参 基因 。 但 TUB 基因 
并 非 适 用 于 所 有 物种 ， 在 甜 标 桃 (Prunus avium) ( 仇 志 浪 等 ，2020) 、 营 草 
(Hemerocallis fulva) Rm, 2020) ~ BÆ (Zantedeschia hybrida) (AWSE, 
20200 不 同 品 种 和 组 织 器 官 中 ，7UB 表达 不 稳定 ， 不 适合 作为 内 参 基 因 。 这 表明 了 内 参 
基因 筛选 研究 应 综合 不 同 物种 、 组 织 器 官 及 不 同 试验 条 件 的 影响 。 
综合 所 述 ， 本 研究 筛选 到 -TUB、p-7UB5 和 a-7TUB3 是 “ 红 元 宝 ' 二 次 花芽 分 化 相关 
基因 表达 分 析 的 稳定 内 参 基因 组 合 ， 能 够 对 开花 关键 基因 进行 较 准 确 的 相对 定量 分 析 ， 该 
结果 为 木兰 属 植物 二 次 或 多 次 开花 的 分 子 机 理 研究 奠定 基础 。 
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